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ty, ° = gre ee 158s 31i29:69 PATUNTy eaeoaTncinacmaR:1 ae 
~TITLE MTHSSINH ; FLoating Point Hyperbolic Sine routine 
IDENT /1-007/ ; File: MTHSINH.MAR Edit: RNH1007 


MARABAAAAAAAAAAALAAAAASALA ALAS AA LEAR EAA EASE EEA RASA LALA AAAS SEERA RAR EADS SS 


COPYRIGHT (c) 1978, 1980, 1982, 1984 B 
ITAL EQUIPMENT CORPORATION, MAYNARD. MASSACHUSETTS. 
RIGHTS RESERVED. 
SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
IN ACCORDANCE WITH THE TERMS OF SUCH LICEN SE AND WITH THE 
USION OF THE ABOVE COPYRIGHT NOTICE. THIS sort ARE OR ANY OTHER 
ES MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
RP NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 


® 
® 
® 
oy 
® 
® 
® 
® 
® 
® 
TRANSFERRED. * 
® 
® 
ww 
® 
® 
® 
® 
® 
® 
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THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
eonPonat io NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


. . . . * oo . . . . * . . . . * 
@aeeneneneneeenenneeeennennannnne 


(RE RRSARAAAAAAARALARLALESELALAALEAASASE EASES ERASERS ERASERS ASRS ALAA AA EAE SS 


fACILITY: MATH LIBRARY 
"ABSTRACT: 
MTH$SINH is a function which returns the floating point hyperbolic sine 


of its single precision floating point argument. The call 1s standard 
call-by-reference. 


VERSION: 01 
ped 

Peter Yuo, 29-Jun-77: Version 01 
MODIFIED BY: 


SOOOOOCOOOOOOOOSOSOOSOSOSOSOSOOSOSOOSOSOOSOOSOOSOOOOOOSOOOOOSoOO 
SOoOOOCOOCOCOSoOOCOOSOOOOSOSOSOOSOOSOSOOSOSCOOSOSOOOOSOOSO 
OOOoooooooooooooooooooooooooooooooooooo 

CON AUS WN $$ OOO NAULS WIN OOO NAUES WN O OONO UE WH OOONOU ES win 
Bete Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Se Ge Ge Ge Ge Se Ge Ge Ge 


SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
PPP RE SS SWNWNIWIWIWININIGIPONOPOPONON NUNN 2 PO Os sO Oe 


ating Po 


A al 
1-00 ; Det 


D 15 
Hyperbolic Sine routine 16-SEP-1984 01:50:41 VAX/VMS Macro V04-00 P 
ed "Torrent Edit History “SEP=1 1386 9 139366 MTHRTL.SRC CIRTHSINH.MAR; 1 _— 


-SBTTL HISTORY ; Detailed Current Edit History 


once 


ALGORITHMIC DIFFERENCES FROM FP-11/C ROUTINE: none 
Edit History for Version 01 of MTHSSINH 
0-2 RTS cel anges. ton to Ms sateeahatateas 


ange gr fondting” mechanism. Put error result in RO before 
tate ng MTHS$SIGNAL in order to allow user modify error result. 


0-3 Six term Taylor series, in powers of the argument, replaced 
by four term vonal ar series, in powers of ARG**2, with 

gxerneng - small arguments; this improves accuracy. 
18-May- $7 Mary Payne 

- Update version number and nBk" 7 EOSvEmA poriee. JBS leony 76 

- ¢ ange MTH__FLOOVEMAT to MTH Ve J03 ee 78 

- Remove SSRADEF macro - not needed. 18s 

- Add "'."’ to the PSECT directive. JBS Deen i 

- Declare externals. SBL 17-May-197 

- han d lower limit for Chebyshev approximation from 2**-11 to 

- Eliminated s second call to EXP for input values between 12.5*|n2 
an 

- Changed all final floating point divisions by 2 to interger 
subrtracts of 1 from the exponent fi 

- Extended maximum range from 87.69 to 128*(n2=88.72. 

poaneged 2° logic for computing EXP(ixi-ln2) to reduce error. 


Change om to G* on call to MTH$$SIGNAL RNH 09-Sept-1981 


SSSI NSAP PAA AA AAAI MMIII 
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Bus’: Beclorar tue Parc of Hednt Gogepatone $4:3700 PATHNTE GaeSericinn mans: 2! cd 


-SBTTL DECLARATIONS 3 Declarative Part of Module 


or 
m 
aon 
cr— 
ro 
De 
> 
ee 


t INCLUDE FILES: MTHJACKET.MAR 


EXTERNAL SYMBOLS: 


GBL 
EXTRN MTHSEXP_R4 
EXTRN MTHS$K_FCOOVEMAT 
EXTRN MTHSSSIGNAL 


: EQUATED SYMBOLS: 


SOoOOoOCoOoOOQOOCOOOCOSOSOOOOOSOOOOSOOOOoOO 


34380 LF_127_LOG_ = *X0F 334380 x 7 ok 
$% 1 LF_128 LOG = maT Def 3 128*ln2e-2%2-24 
420A LF-12.5 L06_2 = “tale +g0A 4 Selng 
04031 LF_LOG_2_HI1 = *X72804031 3: (high 17 bits of Lln2) + 2e*-17 
08700 LF L0G. _LO = *X1100B7D0 : ln2 = LF_LOG_2_HI 
4080 SF_1.0 = *F1.0 3 1.0 
0004 value = 4 3 val 
0 3 


PSECT DECLARATIONS: 


-PSECT _MTHSCODE PIC,SHR,LONG, EXE ,NOWRT 
3; program section for math routines 


SOOoOoooSoSo 


SOOM VGLOOOSOSOOSOSSOSSOSOCOOSSCOSOSOSSOOOOOOOOOOOOOOOSOOOOSOOO 3 


OWN STORAGE: none 


00 
000 CONSTANTS: 
SINHTAB: 
«WORD  *0035120,°0064613 : DECIMAL: 0 1987574E-03 
; a 3$419.2010417¢ : DECIMAL: 83 3320-08 
SWORD *0037452,°012525 > DECIMAL: :1666667E+ 
0 -WORD 0,0 + DECIMAL: 0.00 
SINHLEN = .- SINHTAB/4 


PIRMNPMNNINININ 2 3 FS S$ $$ DDO D000000.0.00.0.0.0.0.0.0.0000 000000000000 


WOdNAUES WIN OOOO NAUE WN OOD NAUES WIN OS OD NA NEW O OO NOAUNE P 


ee ed a ed ed td dd Sd 4 2 2 dd 2s 2 2 2 2 


ve.rf.r 
MACROS: none 
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rout top 1 tt 91:3 741 YAX/VMS Macro V04-00 Page 4 


ng P 
. on Fl 6=SEP-19 1:27:09 CMTHRTL.SRCIMTHSINH.MAR; 1 (4) 


i nt Hyperbolic Sine 
H = Standard Single Precis 


-SBTTL MTHSSINH = Standard Single Precision Floating SINH 


ds + oat 
1-00 


p++ 

2 ; FUNCTIONAL DESCRIPTION: 

SINH = single precision floating point function 
SINH(X) is computed as: 


SIDE EFFECTS: 


Signal: MTH$_FLOOVEMAT if 128*ln2 =< {X! with reserved operand in RO (copied 

to the signal mechanism vector CHFS$L_MCH_RO/R1 by LIBSSIGNAL). Associated 
message is: ‘FLOATING OVERFLOW IN MATH LIBRARY’'. Result is reserved operand -0.0 
unless a user supplied (or any) error handler changes CHFSL_MCH_RO/R1. 


1 
1 
1 
1 
1 
1 3 
1 3 
1 3 
! 3 
1 ; If 3K! _¢2ee-t2, SINK(O) = X. 
1 3 IF Qee-12 =< 1X} < 0.25, SINH(X) = Chebyshey Series 
1 ; If 6.25 =< Ix! < 13:38 n2, SINH(X) = (EXP(X) = EXP(=x))/2 
1 ; If =< IK! < 127#Ln2, SINH(X) = sign(X)sEXPC IX!) /2 
1 ; If 127#ln2 =< 1X! < 128*ln2, then SINH(X) = sign(X)* EXP(!X!-LOG(2)). 
! $ If 128*ln2 =< ixX!, then overflow. 
: CALLING SEQUENCE: 
: SINH.wf.v = MTHS$SINH(x.rf.r) 
1 : INPUT PARAMETERS: 
00000004 801 LONG = 4 3; define Lonquorg multiplier 
00000004 001 x = 1 * LONG 3; Contents of x is the argument 
3 IMPLICIT INPUTS: none 
001 OUTPUT PARAMETERS: 
4 VALUE: floating hyperbolic sine of the argument 
4 IMPLICIT OUTPUTS: none 
80 COMPLETION CODES: none 
3 
$6 


NOTE: This procedure disables floating point underflow, enables integer 
overflow. 


Be Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Se Ge Ge Be Ge Ge 


ENTRY MTHSSINH, “M<IV, R2, R3, R4, RS> 

' 3; standard call-by-reference entry 
disable DV (and FU), enable IV 
flag that this is a jacket procedure in 


GOOD 0909 SINS SDE DDD DDS DDD TTT EBS BS BS SB EE ENNIO 


WR S ODONAUE WR 9 OONAUE WIN 9 ODNAU EWN (OVO ONAUES WH" OO 


MTHSFLAG_JACKET 
MOVAB G*MTHSSJACKET_HND, (FP) 


set handler address to jacket 
handler 


DWOMMNIWMININOOOCOCOCOCOOCOCOCOCOOCOOOOSOOSOCOOCOCOOCCOOOOCOOOCOOOOOOOOOOOOOOO 
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1 
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1 
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1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
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TOO) ATHSSINH “. Standard Single Sine out Fl 6S oats yi 9539359 MIMRIL. Sh RCM MTHSINH.MAR; 1 suet (ay 


case of an error in routine 

If an —— convert signal to user PC 
and resigna 

R5 = X = @value(AP) 


Bete te Ge 


55 04 § MOVF  a@yalue(AP), R5 
0 MOVF R i 
0 ‘g000 BF AA BICw2 4X80 0, R : RO = 
F80 8F 0 8B CMPW o's 3 td ix with fe $2 : 
iA BGEQ ee Qs .10_0 t branch if Xx! § 4 
. a 
4 'y! al 
; KX! < 0.25 " 
3a80 8F 50 8B CPW #*X3A8 : compare !X! with 2**-12 S\ 
1 BGEQ ee. TO_2m1 : branch if {X! >= 2ee-12 7 


=z 
o 
< 
mn 
aD 
uw 
= 

» 
oO 
. 


~~ 2 Um 0=*xX 
04 RET 3 return with result = Argument 


S Qeent2 =< 1K! < 0.25 


Oooo 
Oooo 


o 
7 
oot co — 
Sss 
SOOoooooccooooooooooooo 
WINN AANA IAI AAAI AAAI AIPOPOPONOPONONO PN 2 2 


PEAT SD NINN NNN OWI OOOO PUIOSO OWWVWwO 


ooo 


GEQ_TO_2m12: 
MULF 


3 0.25 =< IKi < 12.5*ln2 


0 
0 
50 50 44 00 0, RO 3 Get ARG**2 for POLYF 
C1 AF 03 50 55 OO POL YF RO, #SINHLEN-1, SINHTAB ; RO = SUM(Ci*X#*(2*7)) 
0 5 ie Coot tic tent is zero 
50 55 44 00 MULF R5, RO 3; MUL by ARG and then 
50 55 40 004 ADDF R5, RO 3 Odd in ARG with overhan 
04 Be RET 3: return with result in R 
4 3 
oe 3 0.25 =< ix! 
4 3 
04 
04 GEQ_TO -o; 25: 
OF 334380 8F 20 51 004 CMPF RO, #LF_127_L0G_2 ; compare Xi with 127*#ln2 
A % 4 BGTR  GTR_THAN_127_L0G_2 : branch if X > 127*Ln2 
4 3 
2 3 0.25 =< IKi =< 127#Ln2 
4 e 
A123420A 8F 50 51 4 CMPF RO, #LF_12.5_ LOG_2 ; Compare :X! to 12.5*ln2. f greeter. 
1E 14 : BGTR one TERA ONLY ; only one call to EXP is mie 
5 
5 


HAAN. AAAI AAI PIPININPINIPININID 3 9 9 MQDOOCOCOOOOCOOOOOOOOOOOOMMmcccac. 


DOONAN WIN SO OOD NIA UNE WIN 9 OD NAUE WIN 9 OD NAME WIN 9 OD NAME W OOO NOUNS 
@ 


—W CCOdOOO. TH MMH NO 


50 55 50 MOVE R5, RO 3; RO = X 
socowee"t 16 JSB MTHSEXP_R4 3 RO = EXP(X) 
0 oD 6 PUSHL RO 3 push EXP(X) on stack 


MTHSSINH ; Floating Point Hyperbolic Sine r9 ies 1 en 984 7:41 YAX/VMS o V04-00 p 
1 007 arnssin NH ., Standard Siete Precis aut Fl 6S 7138 9 £39; 70 UMTHRTL.§ L.S SATHSINH.MAR; 1 mi (dy 
ie 3 4 MNEGF RS, RO 3 R - 
ae i 4 4S 54 188 ATASEXP Ra i = EXP(=X) 
0 4 4 SUBF3 (SPT+ $ = EXP(X) = EXP(-x) 
‘-¢ 4 SUBW Rei 06 0, kor” : RO = (EXP(X) = EXP(-X))/2 
04 7 rk: RET 3 a with result in RO 
7 46 ; 
9 23 3 12.5 =< txt < 127*in2 
se $8? 
76 50 ONE _TERM_ONLY: 
00000000 ' EF 76 138 MTHSEXP_R4 3 RO XPCIX!) 
4 Ni 7C 3 oe i If X S negative, change 
3 3 Be 27 BTR eBay 3 the sign of EXP(iX!) 
0 55 POSITIVE: ; : RO = sign(X) *EXPCIX!) 
50 0080 8F A2 00 28 SUBW #*x0080, RO 3 RO = sign(X) *EXPCiX!)/2 
04 3 8 RET 
089 5 : 
+94 60 ; 127*ln2 =< {Xi 
aoe 
089 65 GTR_THAN ur a4 _L0G_2 
72164381 8F 50 51 0089 64 nO, me _128_L06_2 ; Compare ‘Xi with 128*ln2. If greater, 
20 «18 009 65 BOEO ; overflow is signaled 
009356? : 
999 68 ; 127*Ln2 =< IX! < 128*Ln2 
09 70 ° 
50 72804031 8F 42 $0 71 SUBF #LF_LOG_2_HI, ; RO = {Xi-(Chigh order bits of lLn2) 
00000000 ‘EF ig 099 a JSB MTHSEXP RZ 3; RO = EXPC(iXi=(high order bits of Ln2) 
51 50 11008700 8F 4 7 MULF3 LF TOG” ~2.L0, RO, R1 : 
50 1 4 74 SUBF R1, RO ~ ; RO = EXP CIX! = Ln2) 
33 3 fg a Fg Pe he 
: branc = 
o> Dw 6 MNEGF RO, RO 3: RO = sign(X) * EXP(iX!-LOG(2)) 
04 23 | 10$: RET 3: return with result in RO 


i 128% n2 =< iX!, error 

4 condition value 

RO = ges = reserved operand -0.0 

goes © 21snal mechanism vector 

CHSt onc H_RO/ R1) so error handler 

con modify the result. 
oignes error and use ne? user's PC 
ependent of CALL v 

4 be me - RO pestered tenn CHFSL_MCH_RO/R1 

g 


Sssss 
SHODOOOOS > > > OOO 
—~no = ~_mT 


ERROR: MOVZBL #MTHSK_FLOOVEMAT, -(SP) 
ASHL #15, #T, RO 


00000000'GF 01 


~ 
@ 


CALLS #1, G*MTHSSSIGNAL 
RET 


oS 
rx 


5 
; 
§ 
; 
g 
5 


AAAAAO 
Mrornon——> > > PS 


Ssss 


MTHSSINH 
Symbol table 


LF LOG_2_L0 
MTHSSJACKET_HND 
MTHSSSIGNAL 
MTHSEXP_RG 
MTHSK_FCOOVEMAT 
MTHSSINH 


PSECT name 


. es 
“MTHSCODE 


Initialization 
ommand processing 


Pass 

Symbol table sort 
Pass 2 

Symbol table output 


Psect synopsis output 
Cross-reference output 
Asseabler run totals 


bytes 
There were 


The worki et Limit was 900 
4456 09 ¢ 


Page faults 


pages) of virtua 
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aes es 9133 :09 CMTHRTL.SRCJMTHSINH.MAR; 1 

p 1 
& 1 
Rp 1 

9R 1 

A 


| 
| 
0 


= OF 834 0 
4 
2 i if 
= DO 
= 00000004 
Reeeeeee x 1 
Reeeeaeee x 
Reeeeenee xX 6 
eeeeeeee x 
000010 RG 1 
B88 ig a : 
= 09000004 
94440444 01 
= 00000004 
baa wore enema wena > 
; Psect synopsis ! 
Allocation PSECT No. Attributes 


00000000 < 0.) 00 ¢ 0.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
000000C2 ¢ 194.) 01 ¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 


Pe 


vie 9808:08-68 S008: B 
| 
sol $3:99:00°99  90:08:00:09 


ages. 
nery were used to buffer the intermediate code. 
pages of symbol table space allocated to hold 22 non-local and 1 local symbols. 


source Lines were read in Pass 1, producing 11 object records in Pass 2. 
page of virtual memory was used to define 1 macro. 


(4) 


MTHS 
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vARetd Macro Run Statistics ° ” Ze iE peste aot 91:30:44 m ney d 


cr 
MTHRTL.SRCJMTHSINH.MAR; 1 


tewmeswresececanzreoccomocemaaawmoe $ 


! Macro Library statistics ! 


eee nme eer ee em meme wee} 


Macro Library name Macros defined 
_S255$DUA28:: [SYSLIBJSTARLET.MLB; 2 genic We 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHSINH/OBJ=OBJ$:MTHSINH MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) +MSRCS: 


“=z 


ORPORATION 
OPRIETARY 


AH-BT13A-SE DI 
VAX/VMS V4.0. CO 


PMENT 
AND P 
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